Born-Oppenheimer approximation for open quantum systems within the quantum 

trajectory approach 
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Based on the quantum trajectory approach, we extend the Born-Oppenheimer (BO) approxi- 
mation from closed quantum system to open quantum system, where the open quantum system is 
described by a master equation in Lindblad form. The BO approximation is defined and the validity 
condition is derived. We find that the dissipation in fast variables benefits the BO approximation 
that is different from the dissipation in slow variables. A detailed comparison between this exten- 
sion and our previous approximation (that is based on the effective Hamiltonian approach, see X. 
L. Huang and X. X. Yi, Phys. Rev. A 80, 032108 (2009)) is presented. Several new features and 
advantages are analyzed, which show that the two approximations are complementary to each other. 
Two examples are taken to illustrate our method. 

PACS numbers; 03.65.Yz, 05.30.Ch 



The adiabatic and Born-Oppenheimer (BO) approxi- 
mations are among the oldest approaches in quantum me- 
chanics P, . The adiabatic approximation tells us that 
for a time-dependent system governed by Hamilto- 
nian H{t), if the system is prepared in one of the eigen- 
state \n{t = 0)) of H{t = 0) at i = 0, it will keep in that 
eigenstate \n{t)) of H{t) at arbitrary time t > pro- 
vided the Hamiltonian H{t) is changed slowly enough. 
The Born-Oppenheimer approximation was first given by 
Born and Oppenheimer in 1927 0, which can be formu- 
lated as follows Treating the slow variables as pa- 
rameters, we first solve the fast variables with fixed slow 
variables. Using these solution we obtain an effective 
Hamiltonian for the slow variables. This effective Hamil- 
tonian contains an effective vector potential induced by 
fast variables. Based on this Hamiltonian we can obtain 
a wave function with the slow variables. Thus the total 
wave function can be factorized into a product of two 
wave functions corresponding to the fast and slow vari- 
ables. This method has been widely used in physics and 
quantum chemistry and becomes a fundamental tool in 
these fields 



Due to the unavoidable coupling of quantum systems 
to its environments, most quantum systems are open and 
dissipativefl^. The dynamics of open quantum system 
can be described by the master equation [23, [25j . It is 
then natural to ask: How to extend these approximations 
from closed system to open system? The adiabatic ap- 
proximation has been extended to open systems in differ- 
ent ways, including the Jordan blocks method in Liouville 
space [ia, [l^] , the effective Hamiltonian approach (28l - [30| 
and in weak dissipation limit (3l| . Although the approx- 
imation based on the effective Hamiltonian approach is 
equivalent to the Jordan blocks method[28|], the effective 



Hamiltonian approach has advantages that the extension 
is straightforward and the effective Hamiltonian is easy 
to give. The BO approximation has been extended based 
on effective Hamiltonian approach in Ref. [3^. In this 
paper, we shall extend the BO approximation in a dif- 
ferent way, which is based on the quantum trajectory 
approach. Compared with our previous work, this exten- 
sion exhibits the following new features and advantages. 
(1) We do not need to extent the Hilbert space. This will 
save the CPU (computing) time and memory. (2) All 
eigenstates of the effective Hamiltonian in the method 
are physical states. (3) The eigenstates are easy to ob- 
tain. (4)It is more accurate to treat the jump terms in 
the master equation than that in our previous method. 

Consider a quantum system with two types of vari- 
ables, a slow one X and a fast one Y. Then we can 
divide the total Hamiltonian of the system H into two 
parts 



H = HsiX) + Hf{X,Y), 



(1) 



where Hg{X) only contains the slow variables X. The 
two types of degrees of freedoms couples together through 
Hf{X,Y). We start considering dissipation in the fast 
variables first. The case of decoherence in the slow vari- 
ables will be discussed later. Assuming the dissipation 
is in the Lindblad form, the dynamics for such a system 
can be described by 



d 



(2) 



where the first term on the right hand side represents 
an unitary evolution while the second term denotes the 
dissipation. Here we assume the dissipative term can be 
arranged into the Lindblad form as 
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where Lk = Lk{Y) is the Lindblad operator relevant to 
the fast variables Y. L^pLj^ denotes the jump term. 
Within the frame of quantum trajectory approach, for 
an initial state \(j){to)), one can write the state after an 
infinitesimal time dt as 

pito + dt)= (l-JI^^Pfe) \M{M+^dpk\<l>k){4'km 

\ k / k 

where dpk — {(j){to)\L'lLk\(l){to))dt and the new states are 
defined by 

, , ^ _ {l-tH,ffdt/h)\<j){to)) 



V^l ~ Efc dpk 
LkWo)) 



\Lkmo))\v 



(5) 



with the non-Hermitian effective Hamiltonian defined by 
Hcs = H—^T%Y^^ L^Lk. Under this description, the sys- 
tem will jump into the state with probability dpk, 
and evolve according to non-Hermitian effective Hamil- 
tonian TJofi with probability f — X^fc'^Pfc- This unravel- 
ing is the so-called Monte Carlo wave function method 
|33l435l |. The difficulty here is that the non-Hermitian 
Hamiltonian ffoff contains two types of variables, we will 
solve this problem by applying the BO approximation in 
the no-jump trajectory only. For the non-jump evolution 
\4>o), the time evolution is given by 



in-\^{t))^H,s\^{t)). 



(6) 



Our aim is to solve this equation with the help of BO 
approximation. To this end, we first rewrite as 
iJeff = Hs{X)+H'jiX, Y) with H'jiX, Y) - HjiX, Y) - 

^'^^L^i^Lk- Obviously, Hj{X,Y) is not Hermitian 
that includes all non-Hermitian parts of i?cff- Taking 
the slow variables X as parameters, we can solve the 
eigenstates for H'f{X,Y). We denote its right eigen- 
states by and the corresponding left eigenstates 

by {ipf^{X)\ with complex eigenvalues En{X). These 
eigenstates satisfy the relations (^A^IV'n) — ^mn and 
(V'n IV'n) = 1 for fixed X. We also restrict our discus- 
sion to the non-degenerate case. In order to solve Eq. © , 
we expand the eigenstate of Hcs in terms of It/j^lX)) as 



N 



|<i>)=5]c„|^„(X))|V.^(X,f)), 



(7) 



Define i?„,™(X) = we can rewrite 

Eq.® in a matrix form as 

ino + nv)^ = E^, (9) 

where Hq, H-p and ip are defined by 



where N is the dimension of the fast variables Y and 
c„ {n — 1, 2, 3, N) are the expansion coefficients. Sub- 
stituting Eq.® into the eigenvalue equation Hce\^) = 
after simple calculation, we obtain an equation for 
the wavefunction of the slow variables 



Hi+Ei{X) 

H2+E2{X) 
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Hn.1 
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(10) 

Here we have omitted the same subscripts for simplic- 
ity. Treating Ji-p as perturbation, we can solve Eq.® 
by virtue of the standard time-independent perturbation 
theory. The solution to the zero-order 
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(11) 



can be obtained by the eigenvalue equation 

^nlfn^k ) - ^n,kWn.k /' 



(12) 



where V-n = {Hn{X) + En{Xj) is the zero-order effective 
Hamiltonian for the slow variables. From these zeroth- 
order solutions, one can obtain the higher-order correc- 
tion accordingly. The condition with which we can ne- 
glect the higher-order correction safely is 



H„ 



fl[0]\ 



E. 



[0] 



E, 



m 



< 1, for all fc',n' ^ k,n, (13) 



where {(fi^l^l, \ is the left eigenstate of non-Hermitian 
Hamiltonian Hn- 

Next we consider the dissipation about the slow vari- 
ables. The method used in this case is very similar to the 
discussion given above. In this case, the Lindblad opera- 
tor is replaced by Xk and it is a function of slow variables 
only, i.e., Xk — Xk{X). In the non-jump trajectory, we 
diivide the non-Hermitian Hamiltonian as HeS = H'g(X)+ 
kf{X,Y), where Hi{X) = H,{X) - ^Ek^l^k. The 
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FIG. 1: (Color online) Schematic illustration of a Fabry-Perot 
cavity with an oscillating mirror at the right end. 

method to find tlie eigenstates and eigenvalues for Hf 
is the same as for a closed system. These eigenstates 
are denoted by |^„(X)) with corresponding eigenval- 
ues En{X). We can handle the slow variables in the 
same way. The only difference is Hn,m{X) in Eq.® de- 
fined as Hr,,m{X) = (V'nl (Hs{X) -^Ek^l^k) IV-m)- 
The zero-order effective Hamiltonian H„ for slow vari- 
ables in non-jump trajectory can be similarly obtained. 
Compared with the closed system, this Hamiltonian in- 
cludes a non-Hermitian (usually anti-Hermitian) correc- 
tion {ipn \ ~ f X]/c -^k-^kli^n), which comcs from the dissi- 
pation. 

According to the discussion given above, we can solve 
the dynamics governed by the non-Hcrmitian Hamil- 
tonian Hcff via expanding the total state as |$) = 

X^n k Cn,fe|¥'|^fc)IV'i^^)- Then the Monte Carlo simulation 
for Eq.([5]) is the following. We divide the total evolution 
time T into several steps. The interval of each step is dt. 
In each step, a random number e which distributes in the 
unit interval [0,1] uniformly is chosen to determine the 
jump process. If e < J2k '^Pk^ the total state jumps into 
the corresponding state according to the corresponding 
Lindblad operator, i.e., for e < dpi, it jumps to \4>i}, for 
dpi < e < dpi + dp2, it jumps to \4>2) ■ ■■■ If £ > J2k ^Pk^ 
it has no-jump process. The system evolves according to 
the non-Hermitian Hamiltonian Hcs and the BO approx- 
imation is used. This process is repeated as many time as 
^step = T/dt, and this single evolution gives a quantum 
trajectory. We can recover the final state of the system 
by averaging over different quantum trajectories. 

In the following, we shall present two examples to illus- 
trate our method. After these two examples, we will give 
a detailed comparison with our previous work First, 
we consider a Fabry-Perot(FP) cavity with an oscillating 
mirror at one end, acting as a quantum-mechanical har- 
monic oscillator. Such system can be described by the 
following Hamiltonian 

H — huja) a — hxa^ ax -I- V —mCl^x^, (14) 

where w is the frequency of the cavity field with the cre- 
ation and annihilation operator and a, respectively. 
m,^,x and p denote the mass, frequency, displacement, 
and momentum of the oscillating mirror, respectively. 



X = w/L is the coupling constant between the cavity field 
and the mirror. L denotes the length of the cavity. Tak- 
ing the cavity dissipation into account, the Lindblad op- 
erator in this example is Li — ^/7a. The non-Hermitian 
effective Hamiltonian for the non-jump trajectory can be 
written as 

iJoff = - ^7)a^a - hga^a{b + b^) + nn{b''b + i)(15) 

where b = /^{x + ^) and g = Usually, 
the characteristic frequency of cavity field can reach the 
order of about 10^* Hz, which is much higher than the 
nano-mechanical resonator frequency 10^ Hz achieved by 
current experiments 12]. Under this condition, we can 
divide this Hamiltonian into two parts as Hcs — Hg + 
Hf with Hs = hVl{b'^b + \) and = fi{ijj - \i-^)a^a - 
figa^ a{b + h^). The eigenstate for the fast variables Hf is 
IV'n ) — Wa): whcrc \na) is the Fock state for the mode 
a, and corresponding left eigenstate (i/'n„| = («a| and 
eigenvalue i?„^ = h{uj — i^)na — hgna{b + b^). Putting 
these into Ea. (fl2|) and following the BO approximation 
process, we obtain the Hamiltonian for the slow variables 
as 

nn,=hn{b''b+^)-hgnaib+b'')+h{uj - ^ij)na. (16) 

This Hamiltonian can be solved by a displacement of the 
Fock state [T^l as 

I'PZ.n,) ^ D{a{na))\nb}, 

where D{a) — e^^"~^" is the displacement operator 
with A — b — a, and a{na) = \ni,) is the Fock state 
for mode b. Note that in this model, the off-diagonal 
elements of the perturbation "H-p is zero, so the BO so- 
lution \^na nt) non-jump trajectory is an exact 
solution. We study the dynamics for such a Hamilto- 
nian according to the method given above and compare 
this solution to the solution obtained by the Runge-Kutta 
method in FigH We choose |$(0)) = i(|0)-|-|l))(|0)-H|l)) 
as the initial state. To make the effects of dissipation 
more strikingly, we choose the parameters as a; = 100J7, 
g — Q.lil in the simulation. In Figl2ja) we study the en- 
tanglement between the vibration of the mirror and the 
cavity field. We choose negativity [sl] as the measure 
of entanglement for mixed state. In the simulation, the 
density matrix is calculated by averaging over different 
runs, i.e., from the state vectors {tpiit)) for the different 
trajectories, the density matrix can be constructed as 

P(^) — Ji^f IV'i(0)(V'i(OI- Then we can calculate the 
negativity for p{t). We find that the Hamiltonian Eg. p^ 
can produce entanglement. The dissipation decreases the 
entanglement gradually, and the larger the dissipation is, 
the faster the entanglement decays. In Figl^b) we plot 
the average value of the coordinate x for the oscillating 
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FIG. 2: (Color online) (a) The entanglement measured by the 
negativity as a function of Q,t. (b) The average value of the 

coordinate {x) (in units of Vj^j^) as a function of Qt. (c) 
and (d) The fidelity between the QTA solution and numerical 
simulation (Runge-Kutta method). Other parameters in the 
simulation are cj — lOOQ, g = O.IQ. The results are obtained 
by averaging over Af = 150 runs. 



mirror as a function of fii. From the figure we find that 
when the system is a closed system, the average value of 
the coordinate for the mirror oscillates with time. The 
dissipation moves the curve left. Similarly, the strength 
of the dissipation determines the displacement. In our 
simulation, we average our results over Af — 150 runs. 
To check the validity of our method, we compare our 
results with the results from the Runge-Kutta method. 
We use the fidelity jsj as the measure of difference be- 
tween two density matrix. For mixed state, the fidelity 
is defined as F{pi,p2) = Tr y/ ^/plpi ^fpi- This fidelity 
reaches 1 when the two states are same. In FigslSJc) 
and (d) we plot the fidelity between the BO solution and 
the Runge-Kutta method as a function of fJi for differ- 
ent 7. It is obvious that the fidelities are always larger 
than 99.9% in our simulation for M = 150 trajectories 
[3^. The error is smaller than 0.1%. This confirms that 
our method can reproduce dissipation dynamics for open 
system efficiently. 

Next we briefly discuss the dissipation in the slow vari- 
ables for this model. In this case, we also assume the 
dissipation is in the Lindblad form, then the Lindblad 
operator reads X\ — y/nb and the non-Hermitian effec- 
tive Hamiltonian for the non-jump trajectory is 
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FIG. 3: (Color online) The same as Fig. [2] for dissipation in 
slow variables. 



nian for the slow variables as 

Hn^ = nno{b''b + ^) - figna{b + b'') + huJUa + ^ihiim) 

with hilo = h^l — ^ihn. This Hamiltonian can be solved 
in a similar way with the displacement ao = The 
calculations are similar to the process where the dissi- 
pation in fast variables is taken into account. The nu- 
merical results for this case is shown in Fig. [3] Two 
differences can be seen from the figure: (1) When the 
dissipation in slow variables is taken into account, the 
entanglement never disappears. (2) The amplitude for 
the average value of the coordinate (x) decrease strik- 
ingly due to the dissipation. We should note that in this 
model, the perturbation T-L-p is zero in both slow and fast 
variables' dissipation case. In general, this condition can 
not be satisfied (see the second example and Ea. (l23)) ). 
then we need to restrict the dissipation in slow variables 
to be weak, because strong dissipation in slow variables 
enlarge the perturbation T-L-p, which leads the BO con- 
dition to be broken down. This can be understood as 
follows, large dissipation rate may accelerate the change 
of slow variables, so that it is hard to distinguish which 
are the fast variables. 

In the second example, we consider a neutron moving 
in a static helical magnetic field. 



2ttz 



2ttz 



B = B{z) = B [sai9 cos , sin 6* sin 



cos ( 



1. 
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iJeff = hua^a - hga^a{b + 5^) -f fiVL{b'^b + -) ~ -ifinb'^ }{n) The Hamihonian for such a system is 

^ • ^lB■3 = HK + Hs, 



Divide it into two parts as iJoff = Hs + Hj with Hs = 
hn{b^b + 5) - ^itiKb^b and iJ/ = huja^'a - hga'^a{b + b^). 
The eigenstates for the fast variables are \tjj) = \na) with 
eigenvalues i?„^ = huna — hgna{b + 6^). With these 
knowledge, we obtain the zero-order effective Hamilto- 



H = H{z) = 



(19) 



(20) 



where a — (cr^, ct", cr^) is the Pauli operators. Taking 
the spin relaxation into account, the Lindblad operator 
is L2 = \fn,a~ . Treat the coordinate as parameters, the 
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and corresponding eigenvalues (in units of fiB) 



FIG. 4: (Color online) The polarization of the neutron along 
the 2 axis as a function of time t (in units of nh/fiB) for 
different dimensionless dissipation rate g. We have set 6 — 
7r/4 and initially the spin is in a state | + |). 
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FIG. 5: Validity measure T{g) as a function of the dimen- 
sionless dissipation rate g. The results have been normalized 
in units of T{g = 0). Other parameters in the figure are set 
to satisfy h/^BM^L = 10"^ and hk:,/fj,BML = 2 x 10"". 



non-Hermitian Hamiltonian for the non-jump trajectory 
can be written as 

1 



cos9- \ig sin 6'e-2'^^*/-^ 
sm Ue'^ 



2'Kzi I L 



cos 6 



, (21) 



where g = nh/fiB is dimensionless dissipation rate. For 
each fixed z, this non-Hermitian Hamiltonian has two 
right eigenstates 



• a 27rzi/L 



.cos-\l) -l-sin-e 



cos^e^-^'/^ 



10)^ 
|0)), 



IV^^) = ^(sin||l> 
two left eigenstates 

(^il=iv(cos|(l|-|-sin|e-2--/^ 
(V'^| = 7v(sin|(l|-cos|e-2-/^(0|) 



E± = ^^16- 5^ -8^5 cos 61 

In the above expressions, the angle a is defined as 

4 sin 6* 



tan a 



4 cos 6* — ig 
and the normalized coefficient is 
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Note that for a non-zero dimensionless dissipation rate 
(7, a is a complex number. In this case, the relation 
sin^ ^ -I- cos^ § — 1 holds while | sin ^ | -I- | cos ^ | = 
1 does not. Put these eigenstates and eigenvalues into 
Ea. ([T2)) . we obtain the zero-order Hamiltonian for the 
spatial variables as 

^2 

Hn = -^(V-zXO'+£^n, 

i'„ = i(^^|V|V^«), n = +,-. (22) 

With these knowledge, we study the population transfer 
among the internal state for the quantum system. Sup- 
pose that we prepare the spin of the neutron in the state 
initially and manipulate the particle moving from 
z = to z = L in a, fixed time interval T = 3 (in units of 
Trh/fiB). Setting 9 = j, we study the polarization of the 
neutron along z axis versus time t with different dimen- 
sionless dissipation rate g, the results are shown in FigH) 
In the simulation we take J\f = 400 trajectories. Some 
feature can be seen from the figure: When the dissipation 
is absent, the polarization along z axis oscillates between 
and 1 as a cosine function of time. The dissipation 
leads the polarization damping to —1 oscillatingly. The 
stronger the dissipation is, the faster of the damping is. 
To measure the validity condition, we define a function 



r(5) = 



E, 



[0] 

n' .k' 



E, 



[0] 



(23) 



where 0„,,„ = -^(2(V.^,|V|^«)V + (V-^, |V2|^«)), to 
characterize the violation of the BO condition. From 
Figl5]we can see that the spin relaxation benefits the ap- 
proximation. This is same as that in our previous work 
pslfs^ . and it can also be understood as that the dissipa- 
tion in fast variables benefit the approximation, because 
it accelerates the moving of fast variables and the differ- 
ence between the two types of variables becomes more 
evident. 

It is time to give a detailed comparison between this 
method and our previous approach [s^. We note that 
the differences come from the two methods itself. This 
leads to the following distinct features: (1) In Ref. [s^ . 
the extension is done by effective Hamiltonian approach. 
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which requires to extend the Hilbert space. In the present 
paper, the extension is done according to the quantum 
trajectory approach that does not require to extend the 
Hilbert space. In addition, in the present paper, if the 
initial state is pure, the state will always pure in the evo- 
lution. (2) In the effective Hamiltonian approach, the ex- 
tension is simple and straightforward, however, the eigen- 
states of the Hamiltonian including the fast variables may 
not be physical states, although it gives a correct dynam- 
ics. For example, the last three eigenstates of Hj in Ref. 
[3^ are not physical states. In quantum trajectory ap- 
proach, the eigenstates are all physical states. (3) The 
complexity is different. The analytical solution for the 
non-jump trajectory is easier than the effective Hamil- 
tonian solution. For exarnple, the last three eigenstates 
for the fast variables Hj [32| is a cubic equation, whose 
solution is complicated. For a high dimensional open 
system, the problem become more complicated. But it is 
relatively easy to solve the eigenstates in this paper, i.e., 
the non-Hermitian Hamiltonian in non-jump trajectory 
can be solved more easily than the effective Hamiltonian 
in Ref. [33 [. (4) The method in present paper is more 
accurate to treat the jump term in the master equation. 
This can be understood as follows. When the dissipa- 
tion in slow variables is considered, the non-diagonal ele- 
ments of perturbation 'H-p can be divided into two parts 

asHn.rn = {ip7i\Hs\tpm)-^'iJ2k{'^ri\XlXk\ipm)- The first 
part is the same as that in closed system while the sec- 
ond part only contains a term of dissipation. Other part 
of dissipation is recovered via the jump process. Com- 
pared with effective Hamiltonian solution, in which the 
perturbation includes all parts of dissipation, the quan- 



tum trajectory solution is more exact. Before closing 
this paper, we emphasize that all the discussions in both 
two methods are restricted to the non-degenerate energy 
level, i.e., we assume the closed system Hamiltonian is 
non-degenerate. However, when the dissipation is taken 
into account, new degeneracy is introduced [ssl j in both 
methods. For example, in the second example in this 
paper, even the original system is non-degenerate, the 
non-Hermitian Hamiltonian Hpff can be degenerate at 
— ^ and g = 4,. In Ref. [32[ the degenerate occurs 
a,t 6 — ^, z — and 5 = 8. Obviously, the degen- 
erate points in both methods are different, so these two 
methods are complementary in this sense, in other words, 
when one method is not available due to the degeneracy, 
we can choose another method. 



In summary, we have extended the BO approximation 
from closed system to open system by the quantum tra- 
jectory approach. An assumption that the dissipation is 
in Lindblad form is required. The BO approximation is 
used in the non-jump trajectory and the dynamics can be 
recovered by the Monte Carlo wave function method. As 
illustrations, we give two examples to detail our method. 
The results show that our method can reproduce the dis- 
sipation dynamics for such systems efficiently. A detailed 
comparison with our previous work is also given and dis- 
cussed. 

This work is supported by NSF of China under Grant 
Nos. 10775023, 10935010 and 10905007. 
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